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Abstract 

We show how generaUzed parton distributions (GPDs) can be determined in the case where 
hadrons are described in terms of their partonic degrees of freedom through solutions of dynam- 
ical equations. We demonstrate our approach on the example of two-quark bound states de- 
scribed by the Bethe-Salpeter equation, and three-quark bound states described by three- and 
four-dimensional Faddeev-like equations. Within the model of strong interactions defined by the 
dynamical equations, all possible mechanisms contributing to the GPDs are taken into account, 
and all GPD sum rules are satisfied automatically. The formulation is general and can be applied 
to determine generalized quark distributions, generalized gluon distributions, transition GPDs, 
nucleon distributions in nuclei, etc. Our approach is based on the gauging of equations method. 
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I. INTRODUCTION 



The structure of hadrons is comprehensively described through their Generahzed Parton 
Distributions (GPDs) Q i, S 3 (for recent reviews see Refs. 0, 13). GPDs are usually 
defined in terms of light front correlation functions {P'\qi3{0)qo,{y)\P) , where the bilocal 
operator q/3{0)qa{y) is a product of quark, gluon, or other parton fields, and \P) is a hadronic 
state with momentum P. However, perhaps the most generally applicable way to specify 
GPDs is through the density matrix defined as |0] 

Pap{P\P,k) = J dSe''-y{P'\Tq^{Q)qa{y)\P)c (1) 

where 'T' stands for usual time (?/°) ordering, and subscript 'C indicates that only connected 
contributions are retained. Eq. (Q) describes the virtual (off-shell) scattering of partons off 
the hadronic state, as illustrated in Fig. ^ Then, for the study of deeply virtual Compton 
scattering (DVCS) where one is restricted to the light front {y'^ = 0), Eq. (P) can be used 
to define the light front distribution function 

Pap{P',P,k) = J ^Pap{P',P,k) = I d'ye^'-y{P'\Tq0{Q)q^{y)\P)c5{y^). (2) 

where k = (fc+,k"'"). A further integration over d^k^ then leads to the usual definition of 
GPDs.i 

On the other hand, for studies where rotational invariance is important, like investigations 
of the shape of a nucleon BHEH, Eq. (P) can be used to define a rotationally invariant 
distribution function on the = surface 

p,^(P',P,k) ^ j ^p^^{P',P,k) = I d'ye^'y{P'\Tqp{0)q^{y)\P)cS{y'). (3) 

Also, the same expression for as Eq. ((H), but where g/3 and qa are nucleon fields and \P) 
is the state vector corresponding to an atomic nucleus, can be similarly used to describe the 
nucleonic structure of nuclei. 

Using Eq. (P as the basic starting point for the description of GPDs and other quantities, 
this paper addresses the question of how to extract the density matrix pa/s in the case where 
the hadron's structure is modeled by a dynamical equation describing the mutual scattering 




FIG. 1: Virtual quark-hadron scattering as described by the density matrix Pq,/3 of Eq. 



Here one needs to replace the time-ordered product of fields by an ordinary product, a step justified for 
the diagonal case P = P' in Ref. , and for the general case in Ref . jsj . 
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of its constituents. As a concrete example, we consider the case of a meson or diquark 
modeled by the Bethe-Salpeter equation, as well as the case of a baryon, modeled as the 
bound state of three quarks, whose bound state wave function is found by solving a Faddeev- 
like three-body scattering equation. For this purpose we consider both a four- dimensional 
(4D) formulation of the three-body scattering equations, and a covariant three-dimensional 
(3D) one using the "spectator approach" of Gross jl^ . 

An important requirement in any extraction of is the preservation of sum rules which 
relate GPDs to electromagnetic form factors. In this regard, the density matrix defined by 
Eq. (dl) must satisfy the sum rule of Eq. which relates pap to the electromagnetic vertex 
function F^. Charge and current conservation then lead to two further sum rules: Eq. ()28p . 
and Eq. ()29|) - we shall refer to these three sum rules collectively as 'GPD sum rules'. In any 
given model, of course, these sum rules are not guaranteed. Indeed, the nonperturbative na- 
ture of the model considered here (nonperturbative solutions of scattering equations), makes 
the task of extracting p^,^, while preserving the GPD sum rules, particularly challenging. 

Some years ago a similar problem presented itself: how to determine the electromagnetic 
currents of hadronic systems described by nonperturbative solutions of scattering equations. 
Here, of course, the currents had to obey charge and current conservation. The solution to 
this problem came with the development of the 'gauging of equations method' jl^, 3 IE 3, 
17, 3]- This method not only solves the problem of extracting charge and current conserving 



electromagnetic currents, but it does so in accordance with strict adherence to theory: the 
external electromagnetic field is attached to all possible places in the nonperturbative strong 
interaction processes defined by the dynamic equation model. A further feature of the 
method is that it automatically takes care of all the overcounting problems that plague 4D 
approaches lE IE lEl- Indeed, the gauging of equations method has since been instrumental 
in the construction of current conserving electromagnetic interactions li^, HE ISl j and 
in enabling careful analyses of the overcounting problem |;24.. ,25.. 26]. 

In the present paper we exploit the fact that the gauging of equations method can be 
applied not only to matrix elements of the electromagnetic current operator, but to all other 
field operators as well. In particular, by applying this method to the case of matrix elements 
of bilocal operators, we are able to construct the density matrix p^p, and thus GPDs, of 
hadronic systems described by dynamical equations. Moreover, just like 'gauging equations 
with a local electromagnetic operator' results in currents where the electromagnetic field 
is attached to all possible places in the strong interaction model, 'gauging equations with 
a bilocal operator' results in GPDs where the two external legs g/j and qa (bilocal field) 
originate from cutting all possible bare quark propagators in the strong interaction model. 
It is this completeness of 'cutting bare propagators' that makes the resulting pap satisfy the 
GPD sum rules. 

More specifically, by applying the gauging of equations method for bilocal fields to a 
dynamical equation describing a hadronic state, we obtain an expression for the density 
matrix p^,^ that obeys the GPD sum rules whenever the input quantities - the distribution 
functions that result from cutting all bare propagators in (i) the dressed quark propagator, 
and (ii) the quark-quark potential - themselves satisfy corresponding sum rules. Because 
the gauging of equations method for bilocal fields cuts all bare propagators in the model, 
the resulting density matrix not only satisfies the GPD sum rules, but it does so in the 
theoretically correct way. 

For the concrete model of a baryon described by Faddeev-like equations, the gauging of 
equations method leads to Eq. (jlSj) and Eq. (jl^ for the density matrix p^,^ in the 3D and 
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4D formulations, respectively. These expressions satisfy the GPD sum rules whenever the 
inputs satisfy corresponding sum rules, and can be used directly for practical calculations. 



II. GAUGING 

By 'gauging' we shall mean the transformation G — > G^ where G is the n-point Green 
function 

= (0|Tg(xi)...g(x„)|0), (4) 
and G^ is the 'gauged' Green function defined as the corresponding {n + 2)-point function 

...,Xn]X,y) = (0|Tg(xi) . ..q{x)q{y) . . .g(x„)|0)c (5) 

where subscript c indicates that no contributions with a disconnected piece (0|Tg(a;)g(?/)|0) 
are allowed. To be definite, we shall refer to quark parton distributions in hadrons and thus 
take all the q{xi), . . . ,q{xn), q{x) and q{y) to represent quark fields; however, it should be 
understood that our discussion applies equally well, for example, to gluon distributions in 
which case q{x) and q{y) would represent gluon fields, and to nucleon distributions within 
a nucleus, in which case q{xi), . . . ,q{xn), q{x) and q{y) would all represent spinor nucleon 
fields. To make the connection with GPDs more clear, it is also useful to write Eq. (0) with 
two of the spinor components made explicit: 

G'^p{xi,...,Xn;x,y) = {0\Tq{xi) . . . qi3{x)qaiy) ■ ■ ■ qixn)\0)c- (6) 

The definition of the gauged Green function given in Eq. (jS)) can be considered as an exten- 
sion to bilocal fields of the definition used in the case of coupling to a local field, for example 
the electromagnetic field, which involves the (n + l)-point function 

G"^(xi, ...,Xn,z) = (0|Tg(xi) . . . q{z)T''q{z) . . . g(x„)|0) (7) 

where J'^{z) = g(z)f^g(z) is the electromagnetic current operator. Indeed the gauging 
method used in this paper is closely based on the one developed for the (n + l)-point 
function of Eq. Q [iH H H Hi . 



We are interested in the case where the Green function G is modelled nonperturbatively 
as the solution to an integral equation of the form 

G = + GqVG (8) 

where Go is a product of dressed single particle propagators, Gq is the antisymmetrized 
version of Gq, and V is the interaction kernel (in the 3 — > 3 processes of main interest 
here, V consists of all possible 3-particle irreducible diagrams). In the 'gauging of equations 
method', the (n + ll-point Green function G^ is obtained by 'gauging' Eq. (jSJ with a local 
(vector) field as [H [H [13, [l^ 



G^' = G^"^ + G'^VG + GqV^'G + GqVG^. (9) 

Similarly, we obtain the {n + 2)-point Green function G^ by 'gauging' Eq. (jH)) with a bilocal 
(spinor) field as 

G^ = G^^ + G^VG + GoV^G + GqVG^. (10) 
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FIG. 2: Pictorial representation of 'U-gauging'. (a) Bare particle propagator, (b) Bare particle 
propagator cut through the middle, (c) Inner cut legs form a bilocal external field. 



The one-to-one correspondence between these two types of gauging is self-evident, with 
each of the above two equations leading to solutions (for and G^) that are of identical 
form. When necessary to distinguish between these two types of gauging, we shall refer to 
the transformation G — > G^ according to Eq. (fTIHl as 'U-gauging', and the transformation 
G ^ G'^ according to Eq. Q as '/i-gauging'. 

Although this method of gauging is designed specifically for nonperturbative ap- 
proaches, its idea is rooted in perturbation theory: for example, to any diagram V 
contributing to G{xi, . . . ,Xn), there corresponds a sum of diagrams belonging to 
. . . , x„; X, y), each of which can be obtained by replacing a bare propagator 
dQ{u — v) = {Of\Tqf{u)qf{v)\Of) in V (subscript / indicates a free field or state) by 
dQ^^I^{u,v; x,y) = {Qf\Tqf{u)qp{x)qa{ii)qf{v)\Qf)c where subscript c ensures that only the 
contribution corresponding to Fig. |2tc) is retained. Thus Eq. (fTUI) is just the statement that 
G^ is obtained form G by inserting operator qi3{x)qa{y) within every bare propagator of 
every Feynman diagram of the theory. This insertion corresponds pictorially to cutting the 
bare propagator into two, as illustrated in Fig. |21 Unfortunately not all contributions to 
the G'^ defined by Eq. can be obtained by gauging G in this way. However, in many 
problems, including those discussed in the present paper, either only diagrams obtained by 
gauging are of interest, or the missing ones are easily taken into account. Two examples of 
diagrams that contribute to but that cannot be obtained by [/-gauging G are given in 
Fig. El 

For those used to the functional integral approach it might be more transparent to de- 
fine U-gauging by a functional derivative with respect to an external bilocal field B{x, y) 
introduced into the action S'[g, A\ (where A represents gluons or other fields) as S[q, A\ —>■ 
S[q,A,B] = S[q,A] + f d^x d'^y qlx)B{x,y)q{y) This modification amounts to the re- 
placement of the inverse bare quark propagator as 

do^i^-y) d-^{x,y)=d^\x~y)-iB{x,y) (11) 

in the perturbation theory generated by S[q, A], in order to obtain the perturbation theory 
generated by S[q,A,B]. The usual Green function, the vacuum expectation of the time 
ordered product of field operators (excluding vacuum loops), can be written in terms of a 
functional integral as 

Gi{x,},{y,}, {zk}) = (0|T . . . g(x,) . . . qiyj) . . . Aiz^) . . . |0) 

^ / Dq DqDA e'^I^'-^l . . . q{x,) . . . q{y,) . . . Ajz^) . . . 

J DqDqDAe^Sig,"^^ ' ^ ^ 
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FIG. 3: Example of diagrams which cannot be obtained by gauging, (a) One gluon exchange 
contribution. Joining the end points a and /3 gives a tadpole diagram which is zero, (b) A three- 
pion (dashed lines) contribution which, upon joining together the end-points a and /3, gives a 2 ^ 2 
diagram which is forbidden by G-partity. 



When modified by the presence of the external bilocal field this Green function becomes 

/ Dq DqDA e'^^"'^''^^ . . . q{x,) . . . q{yj) . . . A{zk) . . . 



Introducing the generating functionals 



/ DqDqDAe^^i^^'^^ 



Z[B] = j DqDqDAe^^^'i^^^^\ Z = Z[0] = J DqDqDAe'^^^'^\ 



(13) 



(14) 



the gauged Green function can then be defined as 



G {{xi\,{yj},{zk}]x,y) = -i Z 



^ .5Z[B]Gs{{x.],{y,],{zk]) 



6Bix,y) 



■,B=0 



-i Z-'Z[B] 



SGb{{x^}, {yj}, {zk}) 



5Bix,y) 



SGisiixi}, {yj}, {zk}) 



6B{x,y) 



+ d{x,y)G{{xi},{yj},{zk})\^ 



■,B=0 



(15) 



B=0 



The last term in the Eq. (P3|) should be discarded as it contributes only to those disconnected 
terms which are forbiden by the meaning of subscript c. Note that there are still other 
disconnected contributions to Eq. (fT3j) which should be kept. So the proper functional 
definition of gauging is just the functional derivative: 



G^i{xi},{yj},{zk};x,y) 



SGsiixi}, {yj}, {zk}) 



SBix,y) 



(16) 



B=0 



In this sense Eq. (Iiup is just a statement of the product rule for derivatives. Note that Eq. 
(ITT|) for the modified bare propagator is obvious, but can be derived from Eq. (fT^ as well. 
The gauged bare propagator can be derived using 



^(iog(xi, X2) _ . 5B{xi,X2] 
6B{x,y) 5B{x,y) 



-i6{xi - x)6{y - X2) 



and 



-agg + doB -.-^ — U. 



6B 



6B 



6B 



(17) 



(18) 
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One finds that 



dQQ{xi,X2;x,y) 
d^{xi,X2;x,y) 

As expected from tlie discussion above, tlie antisymmetrizing contribution, do{xi —X2)do{y — 
x), does not appear in tlie expression for the gauged free propagator dQ{xi,X2;x,y). It 
is also clear from Eq. (fT^ how the gauging of the bare quark propagator corresponds, 
diagrammatically, to cutting the bare propagator into two pieces, as illustrated in Fig. |21 
Finally, Eq. (fTBj) . together with the product rule for derivatives, enables one to see that 
gauging any complicated diagram corresponds to cutting bare propagators entering this 
diagram in all possible ways. 



. 6dots{xi,X2) 
5B{x,y) 

. SdoB{xi,X2) 



doeixi - x)doB{y - X2) 
= do{xi - x)do{y - X2). 



(19) 



B=0 



III. EXTRACTING GPDS FROM THREE-BODY SCATTERING EQUATIONS 



A. GPD sum rules 



As discussed in the Introduction, all GPDs can be obtained from the density matrix pa/s 
defined in Eq. ([T)). In turn, pa/3 can be found from the {n + 2)-point Green function G^^, 
as given by Eq. by inserting a complete set of states on either side of the operator 
q/3{x)qa{y) and then taking residues at the bound state poles corresponding to the physical 
states \P) and \P'). By writing = GF^^G and recognizing that G ~ i'$p^p/{P'^ - M^) 
in the vicinity of the = pole, one obtains 

p^piP',P,k) = j d'ye^^-y{P'\Tqp{Q)q^{y)\P), = ^>p,T^^pmp (20) 

where \E'p is the bound state wave function corresponding to state \P): 

*p(xi,X2,a;3) = (0|Tg(xi)g(a;2)g(x3)|P), (21) 

and 

^0/3 = G ^G^pG ^ (22) 

is the corresponding bound state vertex function. To find G^^ we shall use the [/-gauging 
method described in the previous section. This is in direct analogy to what was done in 
Refs. jisl . IgL 17, E3 to find the electromagnetic bound state current which is given as 



f{P\ P) = (P'|g(0)f^g(0)|P) = ^p,F^^p (23) 

where 

r'^ = G-^G^G'K (24) 

The close similarity between the definitions of G^ in Eq. (jSJ and G'^ in Eq. ((Tj) is embodied 
in the sum rule 

E / (01 /^"/^(P', P, k)r>;^ = f{P', P), (25) 
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which follows from Eq. ()2U|) and Eq. ()23|) . Furthermore, since the bound state current 
satisfies current conservation, 



and charge conservation, 



{P'-P)^f{P',P) = 0, (26) 



f{P,P) = 2QP>^, (27) 



where Q is the total charge of the three-body bound state, the density matrix p^^ satisfies 
two further sum rules, 

- ^)'^ E / /^-/^(^'' ^' ''Ka = 0, (28) 



and 



5^ / (0^ ^' = 2^^'- ^^^^ 

Note that \P) is an eigenstate of the conserved charge operator Q = f d?xJ^{x) with cor- 
responding eigenvalue Q] that is, Q is a physical quantity corresponding to the conserved 
Noether current J^{x) = q{x)r^q{x). Although we take to be the conserved electromag- 
netic current operator, it is clear that essentially the same expressions will hold for conserved 
isotopic vector currents (CVC), conserved axial currents (CAC), and partially conserved ax- 
ial curren ts ( PCAC) for which the right-hand side (RHS) of Eq. (I26|) and Eq. (j28|) would be 



non-zero 



B. GPDs of two-body bound states 

Before discussing the GPDs of three-body bound states described by Faddeev-like equa- 
tions, it is useful to first demonstrate the main ideas of this approach on the simpler case of 
two-body bound states described by the Bethe-Salpeter (BS) equation. In this respect we 
note that GPDs in Bethe-Salpeter approaches have already received some attention in the 
literature 2^ • 



To describe the scattering of two distinguishable particles (e.g., a quark and an anti- 
quark) one can use the integral equation for the two-body Green function G, Eq. (jH)) (but 
with the antisymmetrization superscript "P" dropped from this and subsequent equations). 
Then /i-gauging Eq. (jH)) gives Eq. (0) which can be solved to give [13] 

= GT^G, (30) 

r^^ = r^rf^ ^ + c/-ir^ + y/^, (si) 

where di is the propagator of particle i, is its electromagnetic vertex function, and 
is the gauged two-body potential. The two-body bound state current can then be found by 
taking left and right residues of G'^ at the bound state poles: 

f = ^{d^d2 + did^ + did2V^did2)(t) (32) 

where if) = did24> defines the two-body bound state vertex function (p in terms of the two- 
body bound state wave function ip, and (if = diT^di is the /x-gauged propagator of particle 
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FIG. 4: The two-body density matrix of Eq. H35|) specifying the GPDs of two-body bound states. 

i. This agrees with the result first derived by Gross and Riska jHJ]- An alternative but 
equivalent approach is to start with the Bethe-Salpeter equation for 0: 

(f) = Vdid2(f). (33) 

This equation can then be gauged to obtain 

(f)'' = T{did2Y(l) + (1 + Tdid2)V''did2(j) (34) 

where T is two-body t matrix. Taking the residue at the bound state pole of T one again 
obtains Eq. ()32|) . As shown in Ref. j3l|, the of Eq. (jH^ obeys charge and current 
conservation if the inputs cif and obey both Ward and Ward-Takahashi identities. 

Replacing /i-gauging by ?7-gauging in the above derivations, one obtains the equation 
corresponding Eq. (jH^ : 

p = (f){d^d2 + did^ + did2V^ did2)<p (35) 

where p is the matrix whose components are pa^. It is clear that the p of Eq. (j35p and the 
of Eq. (jniD will obey the GPD sum rule of Eq. ^ if the input pairs (c/f , ), and {V, V^) 
each obey corresponding sum rules. This aspect is an essential ingredient of our approach 
and will be discussed in more detail below. The other two GPD sum rules, Eq. ()28p and 
Eq. (|^. are then automatically satisfied because of the charge and current conservation 
properties of j^. Eq. (jH^ is illustrated in Fig. HI 



C. GPDs of three-body bound states 

The sum rules for GPDs, derived in subection A above, constitute important constraints 
satisfied by the exact theory of strong interactions. As such, it is desirable that the same 
sum rules be also satisfied by models that seek to approximate the exact theory. In Ref. 
17 . Til l we showed that for nonperturbative strong interaction models described by integral 



equations, the gauging of equations method, when applied to obtain and therefore the 
bound state electromagnetic current j^, results in both charge and current conservation 
being satisfied. We shall now show that the same gauging method, when applied to obtain 
and therefore p^p, results additionally in the GPD sum rules of Eq. (|^. Eq. (j2HI), and 
Eq. (j2ni) being satisfied. 

Here we apply the U-gauging procedure to the case of three identical particles, so that our 
expression for pap can be used directly for calculations of GPDs in the nucleon or nucleon 
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FIG. 5: Contributions to the three-body density matrix of Eq. H40|) specifying the GPD's of three- 
body bound states. 

distributions in ^He. We assume that the interaction kernel V defining G [Eq. (jH))] is given 
as a sum of only two-particle interactions, so that 

3 1 

1=1 

where Vi is the (fully antisymmetric) interaction potential between particles j and k (where 
ijk is a cyclic permutation of 123), di is the dressed propagator of particle i, and the 1/2 
is a factor arising from antisymmetry [13] • As shown in Ref. |1S], the //-gauging of Eq. (0) 
then gives the following expression for C^: 

= GT^G, (37) 



i=i ^ ^ 



(38) 



Taking left and right residues of G'^ then gives the bound state electromagnetic current of 
three identical particles: 



i=l ^ ^ 



(39) 



where \l/p = did2d3^p defines the three-body bound state vertex function $p in terms of 
the three-body bound state wave function \I^p. Replacing /i-gauging by fZ-gauging in the 
above gives the three-body density matrix 



p(P', P,k) = ^Yl ^^■^'^ {dfd-'d^' + ^v^d, - ^v,d^^ d,d, $p. 



(40) 



Again, it is clear that the p of Eq. pn]) and the of Eq. (jH^ will obey the sum rule of Eq. 

if the input pairs {d^,df), and {v'^,v^) each obey corresponding sum rules, and that 
the sum rules of Eq. (j2Hl) and Eq. are also satisfied because the of Eq. satisfies 
both charge and current conservation. The contributions to p, as given by Eq. (j4Up . are 
illustrated in Fig. El It is noteworthy that the last term of Eq. (last term of Fig. EJ, 
comes with a negative sign; as discussed in Ref. this subtraction term is necessary to 
remove overcounted contributions that are present in the first term. 

As in the two-body case discussed above, one can give an alternative derivation of and p 
by gauging the bound state equation for the vertex function. This has the advantage that one 
gauges directly the equation that one actually solves numerically; this way, for example, any 
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FIG. 6: Illustration of Eq. H43|l for the three-body bound state vertex function "I?i in the spectator 
approach. The on mass shell particle is indicated by a cross. 

approximations used in the solution of the equation will appropriately be taken into account 
in the gauged result. Indeed, because of the numerical difficulty of solving 4D integral 
equations, one is often interested to perform a 3D reduction of the original 4D approach. If 
one invokes such a reduction, it turns out that it is better to reduce the dimension of 4D 
integral equations first, and then gauge the resulting 3D bound state equations in order to 
deduce the expressions for j^^ and p, rather than gauge at the 4D level first, and then try 
and reduce to three-dimensions the 4D expressions obtained for j^^ and p. This point will 
be examined in more detail in the next section. For now, we proceed with the gauging of 
the bound state equation. 

As is well known from three-body theory, the bound state vertex function $ can be 
specified by writing it as a sum of Faddeev components, $ = $1 + $2 + '^'s, where 

$1 = -tMPi2^i- (41) 

In Eq. pTj) ti is the two-body t matrix between particles 2 and 3 and P12 is the operator 
that interchanges particles 1 and 2. The two-body potential vi and corresponding t matrix 
ti are defined to be antisymmetric {P23V1 = —vi and P23i^i = —ti) so that full antisymmetry 
of the three-body wave function \1/ = did2d3^ is ensured. 

In order to find the expression for p^^, as defined above, we can follow Ref. and 
directly U-gauge Eq. (jH]). We note however that Eq. PT|) is a 4D integral equation (i.e. 
integrations are over independent four-momenta), while for numerical calculations it is often 
more practical to have a 3D approach. For this reason we shall examine one particular 3D 
formulation, the so-called spectator approach of Gross which lends itself particularly 



d{p)= 7, ^ ^ S{p) = 27iA{p)S+{p^-m^) (42) 



well to our gauging procedure. In the three-body spectator approach, two of the particles are 
restricted to their mass shell by the following replacement of the usual Feynman propagator 
dip): 

iA{p) 

pZ _ ^2 _|_ 

where A{p) is a function that depends on particle dressing (A(p) = ^ + 171 for structureless 
fermions) and 6~^{p'^ — m?) is the positive energy on- mass-shell 5-function. We refer to 5{p) 
as the "on-mass-shell particle propagator". Thus to obtain the 3D bound state equation 
corresponding to Eq. (j4H) within the spectator approach, all one needs is to replace one of 
the propagators in the 4D approach, say ^2, by 82: 

$1 = -ti52d^Pi2^i. (43) 

A graphic representation of this equation is given in Fig. IHl (for visual convenience we depart 
from the convention used in all other figures and draw Fig. IHl and Fig. [7| using the same time 
direction as in algebraic expressions: from right to left). Gauging then proceeds similarly 
for Eq. (PT|) and Eq. (jlHI) . As the spectator equation of Eq. (jlHj) involves two types of single- 
particle propagators, d and 5, it represents a more general case than that of Eq. (|¥T|). For 
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(a) 



(b) 



>-<M> 




FIG. 7: (a) The density matrix p^p in the spectator model as given by Eq. (|44)) . (b) The corre- 
sponding bound state electromagnetic current j'^ derived in Ref. |ld| . 



this reason it will be sufficient to discuss in detail the gauging of just the spectator equation, 
and then recover the corresponding 4D results simply with the replacement 5 — *• ci. 

The /i-gauging of Eq. ()43p was carried out in Ref. in order to derive the electromag- 
netic currents of three-body systems in the spectator approach.^ It is therefore clear that to 
obtain the density matrix p^js of Eq. (plj) . and therefore the GPDs of the three-body bound 
state described by Eq. (jlHjl . all we need to do is to replace the /x's by t/'s in the results of 
Ref. in particular, this means that we need to replace the input currents of Ref. [3] by 
corresponding input t/-gauged quantities. Before discussing these inputs, we shall ffist use 
the results of Ref. [3| to write down the expressions that would be obtained by U-gauging 
Eq. (gSl). 

For three-body bound states the expression for p, the matrix made of components pap, 
is given by 

p(P', P, k) = <lf 'Pl25l52rf3tf 52rf3Pl2$f - ^f'^l {5^d:, + 52^) ^12'^>f , (44) 

which follows directly from the corresponding expression for the bound state current [l6l |: 

j^(P', P) = <irPi25i52(i3t^52ci3Pi2$f - ^r^i (^2^4 + 52^) ^12$f • (45) 

These expressions are illustrated in Fig. [7| Note that the last two terms of Fig. [7|^a) do not 
give the full one-body contribution to p as a further contribution comes from the gauged 
propagators inside . To find , we ffist need to specify the spectator equations for ti. 



ti = vi + -vidid'iti] 



ti=Vi + -ti62d3Vi. 



(46) 



By gauging these equations one can express in terms of the gauged potential as 

tf = ^ti {5^ds + S2dl^) + h^d,^ (^1 + ^S^dst^j (47) 



which corresponds to Eq. (27) of Ref. [16[. Note that Vi is the sum of all possible irreducible 
diagrams for the scattering of two identical particles, therefore P23f i = fiP23 = —Vi. That 



^ We note that two recent works 0, [s^ have obtained the same results as that of Ref. Q . 
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FIG. 8: Contributions to the three-body density matrix of Eq. (|52|) specifying the GPD's of three- 
body bound states in impulse approximation. 



is why we do not need to use the symmetrised propagator |(52C?3 + '^2<53) in Eq. pUj) in order 
to satisfy the Pauh exclusion principle. 

Although Eq. ()44|) may be the most practical equation for numerical calculations, with 
the help of Eq. we can also eliminate in favour of the gauged potential 

(48) 

Comparing Eq. ()48|) with the corresponding expression obtained by gauging the original 4D 
equation, Eq. pij) : 

p(P', P, fc) = <lf'4 (44 + 44"") - ^12) $ f +5'r (^12 - 444t^f 44(A2-^)$f , 

(49) 

reveals the prescription di — > 5i, 4 — ^ ^2, c?^ — >• 5f , ci^ ^ 5^ that one should use to 
obtain the density matrix in the 3D spectator approach, Eq. (pSj) . from the corresponding 
4D expression of Eq. ()49p . The corresponding three and 4D equations for the bound state 
electromagnetic current are 

/(P', P) = <lf '5i (52% + 52^^) Q - Pi2^ $ f + (^Pi2 - 5245i<524(Pi2 - ^)$f , 

(50) 

and 

3^{P\ P) = $f'4 (ci% + 4ci^) Q - P12) $f + $r {Pi2 - 44cii<44(Pi2 - ^)*f , 

(51) 

respectively. In the impulse approximation where the gauged potential is neglected, we 
have that 



Pi^p(P', P, fc) = $f 5i (52^4 + 524"") - ^^12) $f 



(52) 



which is illustrated in Fig. |H| This is the full one-body contribution to the density matrix. 
Because of propagator 5i in this expression, particle 1 is on mass shell (of course to the right 
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of operator P12 this on-mass shell particle becomes particle 2). The first term on the RHS 
of Eq. (jS21) also contains the gauged propagator 5^, and, as is evident from the discussion 
below, particle 2 can be off mass shell either to the left or to the right of the bilocal field 
{qpqa) vertex. Thus to calculate this first term, one needs to know $f and $f where only 
one external particle is on mass shell. These can always be determined from the spectator 
bound state vertex functions where two particles are on mass shell by using Eq. (|43p and its 
conjugate, $1 = —^iPi252d^ti. Choosing the momenta of particles 1 and 2 as independent 
variables, we may write Eq. ()52|) in the explicit numerical form 

Pin.p(P',P,fc) = j ^^$f'(Pi,P2 + A)5(pi)5^(P2 + A,p2,A:)d(p3) (^-^12) *f(Pi,P2) 

+ / ^^$r(Pl,P2)5(pi)5(p2)rf^(p3 + A,P3,fc) <ff(Pl,P2) 

(53) 

where A = k' — k, = P — pi — p2, and the momenta which are on- mass-shell are labelled 
with a bar over the top. 

The exact expressions for p in either the 3D spectator approach, Eq. (j48p . or the full 
4D one, Eq. (j49|) . correspond to cutting (as in Fig. |2)) every one of the infinite number 
of bare particle propagators that exist in the nonperturbative expression for the strong 
interaction Green function G, Eq. (jH}. In the same way the corresponding expressions for 
the electromagnetic current j^, Eq. (jSUj) and Eq. (|^. correspond to attaching a photon to 
all the bare propagators in the nonperturbative theory. It is this completeness of the photon 
attachment that guarantees the charge and current conservation properties of j'^. 

It is now evident that in the spectator approach, the p of Eq. (|48|) and the of Eq. 
will fulfill the GPD sum rules of Eq. ((23), Eq. (j2HI), and Eq. ^ whenever the gauged 
input quantities, d and , themselves satisfy corresponding sum rules. Similarly, in the 4D 
approach, the p of Eq. (jl^ and the j'^ of Eq. ()5ip will fulfill the GPD sum rules whenever 
the gauged input quantities, d'^ and v'^ , satisfy corresponding sum rules. 



D. Gauged propagator 

The gauged one-particle propagators d^ and d^ form one of the basic inputs to the 
equations describing the three-body bound state current j^, Eq. (j45j) . and the density matrix 
p, Eq. (jH)), respectively. In momentum space these gauged propagators are defined as 

d''{p',p) = J d^xd'^ye'^P'-''-P y\0\Tq{x)q{0)t''q{0)q{y)\0) (54) 

and 

d';i^{p',p,k) = J rf^xA^^'^e^^^'-"+'"-"-^\0|Tg(x)g;3(0)g«(z)g(?/)|0)e. (55) 

Writing the spinor indices of q{x) and q{y) explicitly, the subscript c in Eq. ()55|) means that 
the disconnected part {0\Tqi{x)qj{y)\0) {0\Tqa{z)qf^{0)\0) is excluded. On the other hand, 
the remaining disconnected part, (0|rgi(x)g/3(0)|0) {0\Tqa{z)qj{y)\0) , does contribute to the 
U-gauged propagator as 

C/3(P',P, ^)disc = i2ny5\k - p)dM da,{p). (56) 
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Similarly to Eq. (j25j) . the two types of gauged propagator are related by a sum rule: 

5^ / = ^^^^ 

In addition, c?^ satisfies the Ward-Takahashi (WT) identity 

ip'-p),d^ip',p) = ie[dip) - dip')], (58) 

and the Ward identity 

d^{p,p) = -^e^, (59) 

where e is the charge of the particle (for isodoublet particle fields e is a 2 x 2 matrix). We 
thus can specify two further sum rules for as 



and 



Also, writing 



df'ip'^p) = d{p')T^'d{p), (62a) 
rf^^(p',p,fc) = d(p')ry(p), (62b) 

where and F^^ are one-body vertex functions, the WT and Ward identities give 

{p'-p),T^{p',p) = ze[d'\p') - d-\p)] (63) 

and 

r-(j>,p) = ^e^^ (64) 
respectively, while Eq. (j^7j) implies that 

5^/(0 = ^'^^ ^^^^ 

Note that, to save on notation, we use the same symbols to denote the one-body vertex 
functions of Eq. (fU^ and the three-body bound state vertex functions of Eq. and Eq. 
(|24p - the type of vertex function meant should be clear from the context; moreover, it is 
easy to see that Eq. in fact holds true for both cases. Combining Eq. and Eq. 
with Eq. (jU^ gives the sum rules 

^P' - P)/^ E / 7^ ^Up'^P^ = '^Id-\P') - d-\p)l (66) 
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E/(0iC(P.P.^)a = -^. (67) 

which relate the vertex function to the propagator d. 

In the absence of dressing, in which case we write Eqs. flfj2j) as 

<(p',p) = rfo(p')r^c?o(p), (68a) 
do%(p',p,fc) = cio(p')ro%4(p), (68b) 

where do = i/{'^—m) is the bare particle propagator, and Eq. (j68b|l with all spinor indices 
revealed reads 

doMjpip'^P^ k) = Yl doMp')Kka,ip do,ijip), (69) 

k,l 

one finds that 

= f^, (70a) 
T\i^,^,ifs = i2nrSik-p)5,p5i^. (70b) 

As mentioned previously, the density matrix p, as specified by Eq. (PHj) or Eq. ()49|). will 
satisfy the GPD sum rules only if the gauged inputs d^ and , will themselves satisfy 
corresponding sum rules. For the case of d^ , the above discussion shows that, in the exact 
theory, d^ does indeed satisfy such sum rules, namely, Eq. (jFfjl . Eq. (jHUjl . and Eq. (jHH). 
However, what is of practical interest in our approach is to construct strong interaction 
models of d^ such that Eq. (jKTj) . Eq. (jHIUl . and Eq. (jF)T|) are still satisfied. In this respect, 
one might think of modeling d^ by just its disconnected part, as given by Eq. 

<%(p',p,fc) ^ (27r)V(fc-p)ci.Mp')rf«.(p)- (71) 
In this case the left hand side of the sum rule of Eq. (jKTjl would give 

^ y (2^ = dWd{p). (72) 

Since f'^ ^ F'^, the RHS of Eq. (|7^ is not equal to d'^{p' ,p), so the sum rules of Eq. (|F7jl . 
Eq. (jHUIl and Eq. (jHT)) will not be satisfied; indeed, only in the case of no dressing, for which 
Eq. ()70ap holds, does such a model for d^ satisfy the required sum rules. 

This shows that the connected part of d^ needs to be taken into account if we want the 
GPD sum rules for p to be satisfied in a model with dressed propagators. One way that 
this can be achieved is to construct a purely phenomenological d^ which satisfies the sum 
rules of Eq. (jSTj), Eq. (lUUj). and Eq. However, a theoretically more rigorous way would 

be to first construct a model for the dressed propagator d using the Dyson- Schwinger (DS) 
equation, and then U-gauging the DS equation. The details of constructing d^ in this way 
will be discussed elsewhere (although Fig. 121 discussed later, provides a graphical example 
of such an approach) . 
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E. Gauged on-shell propagator 



A further basic input to the equations describing the three-body bound state current j^, 
Eq. (jlSj) . and the density matrix Eg. (HIl). a re the gauged on- mass-shell propagators 
and 6^ . As discussed in Refs. j3l|, [3], and |l6l|. the /i-gauged on-mass-shell propagator 6^, 
defined as 

5^ip',p) = 2mAip')T^ip',p)Aip) ^^^P" ~ "^^^'^J,^^' ~ ""'^ , (73) 

satisfies both the WT and Ward identities: 

{p'-p),6''{p',p)=te[6{p)-6{p% (74) 

and 

S^{p,p) = -^e^, (75) 



repectively (an explicit proof is given in the Appendix of Ref. (l^). In Eq. (ffH|) . A{p) is the 
factor appearing in the dressed propagator of Eq. (jl^ . and r^{p',p) is the electromagnetic 
vertex function which satisfies the WT and Ward identities of Eq. and Eq. It is 

worth noting that in Eq. ()73|) one can write 

5+(p2 -m2)A(p) = 6+{p^ -m^)Z{^ + m) (76) 



where Z is the renormalization constant of the off-shell dressed propagator d \il 

Eq. ()73p shows how 6'^{p',p) should be constructed if one knows the vertex function 
r'^(p',p) which is derived from the usual off-shell propagator through /i-gauging, see Eq. 
In a similar way we can define the t/-gauged on- mass-shell propagator 6^ as 



S^,{p',p, k) = 27r.A(pOr^,(p',P, mp f^P" (77) 



where vertex function T^^^^lp' ,p, k) is derived from the usual off- mass-shell propagator 
through [/-gauging, see Eq. ()62b|) . Then the sum rule for vertex functions, Eq. (I65|) . implies 
the sum rule for gauged on-mass-shell propagators: 

5^/(0 ^"^^^ = ^^^^ 

In turn, the WT and Ward identities for 5'^, Eq. (f7^ and Eq. ()75|). give the remaining sum 
rules 

(P'-P)^I1 I 7^^«/3(p',P,^)a = ^e[5(p)-5(p')], (79) 

and 

Y.J§yAMv,m.--,e'-M. (80) 
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F. Gauged potential 



The last inputs needed to be considered are the gauged two-body potentials v'^ and f ^. 
As for the gauged propagators, their construction will depend on the nature of the model 
chosen for potential v. However, independently of the model chosen for f, it is a requirement 
of our approach that is constructed so that it satisfy the WT identity 

%v^{PiP'2^PiP2) =i[eiv{p\ - q,P2]PiP2) - v{p\p'2]Pi + g,P2)ei 

+ e2f (p'i,P2 - ?;PiP2) - v{p[p2;puP2 + q)e2] (81) 



and the Ward identity 

v^{p', P — p';p, P — p) = — i 



dv{p', P — p';p, P — p) dv{p', P — p';p, P — p) 

^ dp'^ dp^ ^ 

+(e, + e, f"^"'''' -"A (82) 

M J 



where Cj is the charge of particle i. In the case where v is defined through a finite sum of 
Feynman diagrams or through a dynamical equation, applying the //-gauging procedure to 
the equation defining v will give a that automatically satisfies Eq. ()81|) and Eq. (jH21)- 

With an appropriately constructed u^, all that is left is to construct the ?7-gauged po- 
tential so that it satisfy the sum rule 

Zl J J^'"a/3ip'lP2,PlP2,k)T''p^=v''{p[p'2,PlP2)- (83) 

Again, if v is defined through a finite sum of Feynman diagrams or through a dynamical 
equation, applying the L'^-gauging procedure to the equation defining v will give a that 
automatically satisfies Eq. (j83|l . 

In the case that v is not given in terms of Feynman diagrams or a dynamical equation, 
one would need to construct and ad hoc, in order to satisfy Eq. (|HT|). Eq. (|H^. and 
Eq. (jHl. 



IV. DISCUSSION 

A. Convolution formula 

In the impulse approximation, the expression for the two-body density matrix, Eq. (|35|). 
becomes pimp = (I>{d^d2 + (ii(i^)0; written out in full, 

Pi^p(P',P,A;) = J -0^^p,(P-p)d''{p + A,p,k) d{P-p) MP-P) 

+ I 0^Jp'{p + A)d{P-p)d^{p + A,p,k)Mp)- (84) 

One can use this expression to define a parton-hadron scattering amplitude M{P', P, k) by 
replacing d^ {p + A,p, k) with its disconnected part, as in Eq. (|^. and then chopping off 
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the resulting two parton propagator legs d{k) and d{k'), thus: 

M(P', P,k) = I ^ ^AP - p) {2T:f5\k - p) d{P - p) MP - P) 

+ I ^Jp'ip + ^)diP-p) i2ny5\k-p)Mp) 

= 0p,(P - k) d{P - k) (f)p{P + ^p'{k') d{P - k) (j)p{k). (85) 
Eq. ()84|) can then be written as^ 

Pi^p(P', P,k) = J ^ M(P', P, q) + A, q, k) (86) 

which expresses pimp as a convolution of the amplitude M and the one-body distribution 
function d^ . 

In a similar way, the three-body density matrix in impulse approaximation, Eq. (j^Hj) . is 
given by the convolution formula: 

A^p(P', P,k) = J ^ Ms{P', P, q) 5^{q + A,q,k) + j ^ M,(P', P, g) rf^(g + A, q, k) 

(87) 

where 

Ms{P\ P,k) = j ^^^'r(Pi,P2 + A)5(pi)(2vr) V(fc - p2)rf(P3) Q - $f (pi,P2) 
'^'^\<lf'(pi, A;')5(pi)rf(P - Pi - A:) Q - P12) $f (pi, fc) (88) 



(2vr) 



and 



Ma{P\P,k) j ^,^fi{\pi.P2)KPi)KP2){27:Y6\k-p,) <ff(Pl,P2) 



(271)4*^^ (Pi, ^ - Pi - mPi)KP - Pi - A:) " 

(89) 

It is worth emphasizing that in both the two- and three-body convolution formulas above, 
the density matrix pimp{P',P,k) will satisfy the GPD sum rule of Eq. (j25|l . as long as 
the input one-body distribution functions d^ and 6^ satisfy the corresponding sum rules 
of Eq. ^7\i and Eq. (f75|) . respectively. In turn, one can ensure Eq. (jHTj), as well as the 
important constraints of Eq. (jUUI) and Eq. (jUTj) . if d^ is constructed by applying the U- 
gauging procedure to the strong interaction model used for the one-particle propagator d ; 
that is, if d^ results from cutting all the bare propagators existing inside d (note that the 



With spinor labels exposed, Eq. ® reads p'^^iP', P^k) ^ J A , ^»j(-P', P. q) dYa.jBi<l + A, q, k) 
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FIG. 9: (a) The dressed propagator d in rainbow approximation, (b) The distribution function 
dF = dV^ d obtained by cutting ah possible bare propagators in the rainbow perturbation series 
for d shown in (a). The propagator d defined by (a) and the vertex function defined by (b) will 
satisfy the GPD sum rules of Eq. Eq. and Eq. (|H7jl . 

propagator d is used in a convolution formula in three places: (i) in the construction of the 
bound state vertex function 0, (ii) in the integral defining the amplitude M, and (iii) in 
the expressions = dV^d and d^ = dV^d). To illustrate this point, consider the so-called 
rainbow approximation for d illustrated in Fig. El^a). It is easy to see that the d^ obtained 
by cutting all the bare propagators of this d is given as in Fig. Efb). An essential feature of 
the ladder sum of Fig. El^b) is that it is constructed from just the same dressed propagator 
d that was gauged to obtain d^ . Thus, if the d of Fig. [Ufa) and the d^ of Fig. IHl^b) were to 
be used in a convolution formula simultaneously, the GPD sum rule of Eq. would be 
satisfied. If, on the other hand, one were to construct a ladder sum as in Fig. El^b) but with 
hare propagators d^, one could not also use d^ in the convolution formula without destroying 
the GPD sum rule of Eq. ()25|) : indeed, a very special dressed propagator would be needed in 
the convolution expression (one that when [/-gauged yields Fig. Efb) with bare propagators) 
in order to satisfy the sum rule. 

B. Valence quark contribution 

In subsection III C, we derived three equivalent expressions, Eq. fl40|) . Eq. (jl^ and 
Eq. (jlSjl . for the density matrix p^p of a three-body bound state within the 3D spectator 
approach. This pap satisfies the GPD sum rules if the input distributions d^ and satisfy 
corresponding sum rules (note that, by construction, if d^ satisfies the sum rules then so 
does 5^). This result was achieved by first reducing the corresponding 4D equation, Eq. 
(I4ip . to three dimensions, by setting two of the three particles to be on their mass-shell, 
and then gauging the reduced equation. It is therefore interesting to note, that in order to 
determine the valence quark contribution to a bound state current or GPD, what is often 
done is the equivalent of reversing the mentioned ordering: first gauging a 4D dynamical 
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equation, and then doing a 3D reduction by putting particles on-mass-shell. 

To illustrate this, consider the valence quark contribution to the light-front distribution 
function, Eq. (0), for a bound state of two bare quarks. For simplicity let's also use the 
impulse approximation. Using Eq. (jSH), which was obtained by gauging the 4D two-body 
bound state equation, Eq. (dH), and the fact that d^ip'^p^k) = (27r)^(5'^(p — k)do{p')dQ{p), 
we get 

[4>p.{P - k) doik + A) do{k) do{P - k) MP - k) 
+M{k') do{k + A) do{k) do{P - k) 0p(A;)] . (90) 

Although we started with a two-body equation, this expression contains more than just two 
valence quark contributions (unless the two-body potential V is structureless). To single 
out just the valence quark contribution, we make the replacement 

d^{P - k) = d^{p) ^ 5o(p) =27r(^ + m)^(p+)5(/-m2), (91) 

which corresponds to picking up only one singularity when integrating over k~ . The ignored 
singularities correspond to the non-valence quark contributions. 

Similarly, in the three-body case, gauging of a 4D equation leads to Eq. (jlUj) . which in 
impulse approximation gives 

Pi^p(P',P,k) = ^ I ^^^^p,{p[,p'^)d^(plp,,k)Do{p,,p,)^p{puP2) 

i=l 

^^^^^P'{p[,p'2)d^{p:,P.k) [DoVDo]{p'^,p„p,)MpuP2) 
Ztt Ztt Ztt Zn ■' 

(92) 

where Do{pj,pk) = do{pj) do{pk). We note again that the last therm of Eq. (jU^ is a subtrac- 
tion term that removes the overcounted contributions present in the preceding term. The 
valence quark contributions are then found by taking residues at all two-quark poles in the 
(jk) subsystem. It is not difficult to see that sum of these valence contributions corresponds 
to formally replacing the propagators do{pj)do{pk) by corresponding on- mass-shell propaga- 
tors So{pj)6o{pk) in the first term of Eq. and then dropping the subtraction term - the 
resultant expression does not overcount even though it may seem that the subtraction term 
has been neglected. 

Although such expressions for valence quark contributions are popular in the literature, 
they correspond to gauging 4D equations first, and then setting particles to their mass shells, 
and thus they do not satisfy GPD sum rules. 



C. Final comments 



In this paper we have demonstrated our ^/-gauging method by obtaining the density 
matrix Pq,/3 of Eq. (jT)), and therefore GPDs, for the specific cases of a two-quark bound 
state described by the Bethe-Salpeter equation, and a three-quark bound state described 
by covariant 3D and 4D Faddeev-like equations. The method corresponds to cutting all 
possible bare quark propagators in the strong interaction model, and ensures that GPD 
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sum rules are satisfied automatically. Also, because both our 3D and 4D formulations are 
Lorentz covariant, the derived GPDs obey polynomiality However, it is important to 
emphasize that our procedure is general, applying to any type of dynamical equation, as well 
as to various types of distributions. For example, the [/-gauging procedure can be applied to 
effectively cut all the bare propagators existing within a given dynamical equation model, not 
. ust those of quarks. In this way one could derive the model's generalized gluon distribution 
^, pion distribution, etc. With sum rules corresponding to current conservation being 
ensured, one can, for example, selectively turn on and off [/-gauging of different constituents 
in order to study the fraction of a hadron's quantum number (e.g. spin) carried by a quark, 
gluon, meson, etc. 

The whole approach holds also for nucleon distributions in nuclei, e.g., when the spectator 
equation of Eq. (pSj) is for three nucleons. Similarly, one could consider matrix elements as 
in Eq. (Q), but taken between three-particle states other than {P'\ and \P). Apart from 
these three-body bound states, one can have states of three free particles and those where 
one particle is free while the other two form a bound sub-system. All these transitions are 
considered in Ref. [HI [igI [itI . [isl for the case of electromagnetic currents. Just as was done 



for the case of the bound state current of Eq. (jl^ above, the results for transition currents 
m Refs. [HHI llTl Il8l | can be used directly to find GPDs simply through the replacements 



df'^d^, v^^v^, S^^S^. (93) 

Off-diagonal transitions —>■ nN in context of GPDs have also been suggested and 
their importance emphasized in Ref. In this regard we should mention that applying 
techniques of the present paper to the vrA^A^ equations derived in Ref. [l^, , and their 
electromagnetic currents [20|, would lead to a comprehensive description of the GPDs for 
the transitions A^A^ — *■ 7rNN,d —>■ nd, etc. 

With the same techniques, mesonic correction to nucleon GPDs can be calculated. For 
example, within the NJL model of the nucleon, one can use the static approximation where 
the mass of the exchanged quark is taken to infinity. Then t he q uark-diquark Green function, 
which enters the expression for the mesonic corrections |34|, is calculated algebraically, 
making the calculation of corrections very practical. In Ref. j3J|, we have already shown 
how to calculate mesonic corrections to electromagnetic currents of the NJL model, and 
using the [/-gauging method proposed here one could calculate the mesonic corrections to 
GPDs so that all sum rules are satisfied. The important point is that our approach enables 
these mesonic corrections to enter calculations of GPDs and currents consistently. 
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